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ABSTRACT 

Foam  generation  is  a  promising  method  to  increase  the  efficiency  of 
enhanced  oil  recovery  processes.  The  present  study  began  with  the  obser- 
vation that  the  pinch  off  of  bubbles  upstream  of  an  obstacle  (pore  throat) 
is  a  mechanism  for  bubble  (i.e.  foam)  generation.  We  find  support  for  this 
mechanism  through  a  study  of  the  laminar  to  slug  flow  transition  in  an 
idealized  physical  and  geometrical  setting.  Transitions  between  multi- 
phase flow  regimes  are  of  basic  interest.  A  surprising  conclusion  of  our 
study  is  the  reversal  of  fingering  stabiUty  caused  by  a  strong  transverse 
flow  field. 
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1.  Introduction 

This  study  is  motivated  by  the  problem  of  foam  generation,  which  is  a  promising 
method  to  increase  the  efficiency  of  enhanced  oil  recovery  (EOR)  processes.  The  foam 
serves  to  retard  flow  in  high  mobility,  or  conduction  channels,  and  it  thereby  increases 
the  efficiency  of  the  EOR  process  in  contacting  low  permeability  portions  of  the 
petroleum  reservoir.  Our  study  began  with  the  observation  [6]  that  the  pinch  off  of  bub- 
bles upstream  of  an  obstacle  (a  pore  throat  in  the  rock  matrix)  was  a  mechanism  for  bub- 
ble formation.  We  find  support  for  this  mechanism,  which  we  examine  in  simplified  phy- 
sical and  geometrical  settings.  As  we  now  discuss,  the  transitions  between  distinct  multi- 
phase flow  regimes  are  also  of  interest  to  fundamental  science  and  provide  a  funher 
motivation  for  the  present  study. 

We  regard  multiphase  flow  regimes  as  phases,  in  the  sense  that  within  each  regime, 
the  constitutive  law  which  defines  the  multiphase  flow  equations  varies  smoothly,  while 
along  the  boundary  between  the  distinct  regimes,  which  may  occur  abruptiy  in  the  state 
space  of  the  flow  parameters,  such  as  the  void  fraction,  the  equilibrium  constitutive  rela- 
tions may  change  discontinuously. 

From  this  point  of  view,  a  multiphase  flow  problem  we  have  studied  extensively, 
namely  Rayleigh-Taylor  instability  [3,5],  can  be  regarded  as  a  metastable  transition 
between  bubbly  and  slug  flow.  The  study  of  the  bubble  envelope,  or  the  outer  layer 
between  the  mixing  layer  and  the  ambient  heavy  fluid,  is  dominated  by  bubble  merger, 
which  is  the  microscopic  mechanism  for  the  bubbly  flow  to  slug  flow  transition.  The 
effort  to  discuss  metastable  phenomena  using  only  equilibrium  state  variables  leads  typi- 
cally to  complex  characteristics  and  imaginary  sound  speeds,  while  a  direct  analysis  of 
nonequilibrium  thermodynamics  leads  to  new  dependent  variables,  whose  dynamics  is  a 
rate  law.  For  the  Rayleigh-Taylor  problem,  the  new  def>endent  variable  is  the  bubble  size 
at  the  bubble  envelope.  The  formulation  of  a  rate  law  lies  at  the  heart  of  the  proposed 
statistical  models  [9, 10, 14]  and  is  also  closely  related  to  the  much  studied  formula 
h  =  oAgt^  for  the  acceleration  of  the  bubble  envelope  [7]. 

In  the  present  work,  we  make  no  connection  to  the  averaged  multiphase  flow  equa- 
tions and  we  investigate  instead  the  dynamics  at  the  level  of  individual  bubbles  and 
slugs.  Studies  of  this  nature  have  as  their  objective  the  definition  of  constitutive  laws  and 
rate  laws,  i.e.  the  definition  of  equilibrium  and  nonequilibrium  thermodynamics  of 
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multiphase  flow.  Within  this  microscopic  framework,  we  study  in  detail  a  mechanism 
leading  to  a  transition  to  slug  flow.  Our  approach  is  to  seek  a  consistent  picture,  combin- 
ing simple  theoretical  arguments  with  first  principle  computations  based  on  the  front 
tracking  method  [4].  Other  approaches  to  multiphase  flow  are  presented  in  [1, 2, 12, 13]. 

The  simplified  physics  we  consider  is  that  of  miscible  displacement  in  two  dimen- 
sions, i.e.  Hele-Shaw  flow.  The  equations  for  two  fluid  Hele-Shaw  flow  are 

v=-^-^VP,     i=o,s,  (1.1) 

Vv=0. 

Here  the  two  fluids  are  denoted  oil  (o)  and  solvent  (5),  while  |i,  is  the  viscosity  of  the  oil 
or  solvent.  The  boundary  between  the  two  fluids  is  assumed  sharp  and  it  moves  with  the 
velocity  v.  v  is  discontinuous  across  the  interface,  but  its  normal  component  is  continu- 
ous, so  that  the  interface  motion  is  well  defined.  P  denotes  pressure. 

It  is  well  known  that  these  equations  are  unstable  on  all  length  scales  [8]  and  can 
develop  essential  singularities  in  finite  time  [11].  In  this  sense,  microscopic  physics  (sur- 
face tension  or  capillary  pressure)  will  regularize  the  flow.  Often,  and  this  is  the  point  of 
view  assumed  here,  the  regularizing  length  scales  are  small  relative  to  computational 
length  scales,  and  so  regularization  can  be  safely  achieved  on  numerical  grounds.  This 
small  length  scale  limit  is  appropriate  to  flow  on  a  reservoir  or  inter-well  spacing  length 
scale,  but  is  not  the  correct  description  of  flows  on  the  scale  of  laboratory  rock  samples  or 
of  rock  pores.  In  this  sense,  our  equations  are  a  simplification  of  the  relevant  physical 
equations. 

The  geometry  as  well  as  the  physics  is  simplified  in  this  study.  We  consider  in 
order  of  increasing  complexity,  flow  in  a  uniform  channel,  in  a  channel  of  variable  cross 
section,  around  a  simplified  obstacle,  and  around  a  random  array  of  obstacles. 

A  surprising  conclusion  of  this  paper  is  the  reversal  of  fingering  stability  caused  by 
the  strong  transverse  flow  field.  In  other  words,  we  find  fingers  formed  when  the  heavy 
(more  viscous)  fluid,  i.e.  oil,  displaces  the  light  fluid. 


2.  Initiation  of  Instability 

After  simplification,  the  problem  of  flow  in  a  uniform  channel  has  three  dimension- 
less  parameters.  These  are  the  mobility  ratio  M  =  — ,  the  ratio  V  =  —  of  the  in  flow 

velocities  v,  and  the  fractional  oil  channel  width  W  = .   See  Figure  2.1a.  We 

Wo  +  Ws 

assume  M  >  1 .  This  is  no  loss  of  generality,  as  the  opposite  case  is  handled  by  an  inter- 
change of  the  two  fluids.  There  is  no  flow  through  the  upper  or  lower  channel  boun- 
daries, the  in  and  out  flow  rates  at  left  or  right  boundaries  are  assumed  to  be  constant  in 
each  fluid,  and  the  total  in  flow  and  out  flow  must  be  equal  by  incompressibility.  The 
total  flow  is  a  constant,  which  can  be  removed  from  the  problem  by  rescale  of  the  time 
variable.  Thus  the  boundary  conditions  are  set  by  the  ratios  of  the  in  flow  and  out  flow 
fluid  velocities.  This  simple  geometry  is  modified  by  a  small  barrier  near  the  inlet, 
separating  the  two  flow  channels.  The  equation  MV  =  1  defines  an  equilibrium  (laminar) 
solution.  The  in  flow  velocities  or  the  dimensionless  ratio  V  is  taken  to  be  arbitrary,  while 

the  out  flow  velocity  ratio  Vo^i  =  — >  and  hence  the  out  flow  velocities  themselves,  are 

assumed  to  be  determined  from  the  equation  MVout  =  1- 

The  regime  MV  >  1  corresponds  to  excess  oil  flowing  into  the  channel,  as  is  shown 
in  Figure  2.1,  while  the  regime  A/V  <  1  corresponds  to  excess  solvent.  In  either  case,  for 
MV  ;t  1,  the  initial  conditions  do  not  corresf>ond  to  a  solution  of  the  flow  equations,  and 
the  resulting  transients  are  the  object  of  this  study. 

To  understand  the  details  of  the  pinch  off  phenomena,  we  first  consider  the  r  =  0 
velocity   field,   v  (r  =  0).   This   field   can   be  represented   as   a   sum   of  two   terms, 

V  =  v'  -I-  (v  -  v'),  the  first  of  which  is  a  uniform  flow  field  v',  coinciding  with  the  out  flow 
(right)  boundary  conditions.  This  field  is  purely  horizontal,  and  is  M  times  as  large  in  the 
solvent  as  in  the  oil.  The  second  term,  v  -  v',  is  a  nonuniform  flow  field,  which  has  zero 
out  flow  boundary  conditions,  and  has  net  in  flow  of  oil  and  net  out  flow  of  solvent  at  the 
left  boundary,  if  MV  >  1  .  Otherwise,  if  MV  <  1,  the  in  flow  and  out  flow  at  the  left 
boundary  are  reversed.   From  this  picture,  see  Figure  2.1(b),  the  nonuniform  flow  field 

V  -  v'  causes  entrainment  of  the  fluid  with  excess  in  flow  into  the  other  flow  channel.  In 
other  words,  the  r  =  0  velocity  will  give  rise  to  the  initiation  of  a  finger  into  the  channel 
which  is  deficient  in  flow.    In  contrast  to  the  Taylor- Saff man  theory,  the  finger  is 


transverse  to  the  channel  flow  produced  by  v'. 

To  measure  the  strength  of  the  pinch  off  we  introduce  the  ratio  y  of  the  nonuniform 
flow  to  uniform  flow, 


Y<  = 


V,  -  V, 


V; 


J  =o,  S 


(2.1) 


='V     --^    — . 


■4         i         1 


(a) 


(b) 


Figure  2.1.  The  central  channel  is  solvent  located  between  two  layers  of  oil.  (a)  The 
velocity  field  of  the  total  flow  v.  The  injection  rate  on  the  right  boundary  satisfies 
MV  =  1.  The  left  boundary  velocity  ratio  V  is  larger  than  the  equilibrium  value,  i.e. 
AfV  >  1,  which  results  in  excess  in  flow  of  oil.  (b)  The  velocity  field  of  the  nonuniform 
flow  V  -  v' .  There  is  no  flow  across  the  right  boundary.  The  excess  oil  in  flow  fffoduces 
the  transverse  field  which  drives  the  solvent  flow  backward  at  the  inlet.  As  the  transverse 
field  of  the  nonuniform  flow  is  strong,  a  small  part  of  the  solvent  flow  in  the  computation- 
al domain  may  be  driven  back  into  the  inlet  For  this  reason,  two  small  inlet  barriers 
which  separate  the  oil  and  solvent  in  flows  were  installed  to  avoid  this  situation. 


Then  y^  can  be  expressed  by  terms  of  M,  V  and  W  as  follows 

_il-W).(MV-\) 

To  study  the  nonuniform  flow  v  -  v',  we  examine  the  variation  of  5v^  =  ^^  ~  ^'y  along  the 
interface.  For  two  cases,  W  =  .5  and  in  the  asymptotic  limit  of  large  M,  it  is  possible  to 
determine  the  behavior  of  v^  analytically,  as  we  show  in  Appendix  A,  namely 


v^  =  a:  In 
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\+expi-TVcl2WH) 
l-e\pi-Kx/2WH) 


(2.3) 


Ui  5v  1  —  iio^v  0 
where  K  =  -— — ■ — — — =-  with  5v,  =  v,  -  v/  for  W  =  0.5andK  =  Yov'o  for  large  M. 

In  the  cases  when  (2.3)  is  valid,  we  conclude  that  for  small  x,  Vy  is  approximately 

proportional  to  In  x,  except  that  the  singularity  at  a:  =  0  has  been  removed  by  the  insertion 

nx 


of  the  inlet  barrier  between  the  two  inlet  regions.  For  large  x,  Vy  decays  like  e    ^^^  ,  so 
that  the  nonuniform  vertical  flow  due  to  excess  oil  flow  is  mainly  limited  to  the  region 

0<x  < .  Beyond  this  region,  the  flow  is  approximately  equal  to  its  uniform  value 

2WH 
v'  given  by  the  right  boundary  conditions.  Therefore  the  length  scale measures  the 

of  horizontal  penetration  distance  of  the  nonuniform  flow.  In  Figure  2.2,  we  plot  v^  vs.  x 
along  the  interface,  and  display  the  effects  of  change  in  M,  W  and  V  on  v^. 


3.  Instability  Growth 

Fingers  growing  into  a  strong  transverse  flow  have  qualitative  features  not  present 
in  the  classical  Taylor-Saffman  instability.  In  this  section,  we  construct  a  simple  model 
of  a  square  wave  finger  in  a  strong  transverse  flow  field  and  give  an  approximate  analysis 
to  predict  the  main  features  of  this  flow.  We  then  present  a  study  of  the  exact  solution  of 
the  Hele-Shaw  flow,  for  fingering  into  a  transverse  flow  field,  and  use  the  simple  model 
analysis  to  understand  the  observed  behavior  of  the  computed  solutions  of  the  full  Hele- 
Shaw  equations. 

The  initiation  of  fingering  was  discussed  in  the  previous  section,  and  is  caused  by 
the  nonequilibrium  (A/Vu,  *  1)  inflow  boundary  conditions.  We  distinguish  two  cases. 
For  excess  inflow  of  the  more  viscous  fluid  (oil),  we  find  a  series  of  narrow  fingers, 
repeating  in  a  quasi-periodic  or  chaotic  pattern  for  an  extended  time  period.  Depending 
on  their  length  and  the  width  of  the  other  channel,  the  fingers  may  or  may  not  grow  to  an 
amplitude  which  causes  pinch  off,  i.e.  slug.flow,  as  is  discussed  in  §4.  For  excess  inflow 
of  the  less  viscous  fluid  (solvent),  we  find  a  smooth  transition  to  a  laminar  flow  in  which 
the  solvent  channel  widens  as  one  moves  downstream,  and  approaches  its  equilibrium 
width,  i.e.  a  width  for  which  the  resulting  fluid  velocities  satisfy  MV  =  1. 
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The  patterns  in  these  two  cases  are  somewhat  surprising,  since  narrow  fingers  are 
produced  in  the  Taylor- Saff man  stable  case,  while  wide  fingers  or  nonfingering  behavior 
occurs  in  the  Taylor- Saff  man  unstable  case. 


(a) 


(b) 


(c) 


Figure  2.2.  The  velocity  component  v^  vs.  x  plotted  along  the  r  =  0  interface  for  the 
internal  0<x/H  <  1.5.  Each  plot  displays  the  effect  of  variation  of  one  of  the  three 
parameters  M ,  W  and  V .  The  two  curves  are  rescaled  in  each  case  to  have  the  same  max- 
imum height,  (a)  i)  M  =  40,  ii)  M  =  1.  For  both  cases  V  =  1.25,  W  =  0.5.  (b)  i)  W  =  0.9, 
ii)  W  =  0.1.  For  both  cases  V  =  1.25,  M  =  40.  (c)  i)  V  =  0.25,  u)  V  =  1.25.  Fot  both  cases 
M  =  40,  W  =  0.5. 

The  simple  model  replaces  the  finger  by  a  square  wave.  The  dynamics  of  the  square 
wave  is  then  replaced  by  a  one  dimensional  noninteracting  layer  ansatz,  in  which  three 
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horizontal  layers  are  each  treated  as  independent.  The  layers  with  no  finger  have  a  purely 
linear  pressure.  In  the  layer  with  excess  inflow,  the  pressure  will  have  a  slope  with  a 
larger  magnitude.  In  the  layer  with  the  finger,  there  arc  two  cases  to  consider.  For  a  more 
viscous  finger  (oil)  in  a  less  viscous  ambient,  the  finger  is  a  region  of  very  steep  pressure 
gradient.  For  a  less  viscous  finger  in  a  more  viscous  ambient,  the  finger  is  a  region  of 
very  small  pressure  gradient.  In  terms  of  the  horizontal  motion  in  the  layer  and  Darcy's 
law,  the  oil  finger  is  entrained  in  a  rapid  flow,  and  is  moving  rapidly,  so  has  a  large  velo- 
city and  large  pressure  gradient.  The  solvent  finger  is  entrained  in  a  region  of  relatively 
slow  flow,  so  has  small  pressure  drop  across  it  The  results  of  these  pressure  fields  are 
shown  schematically  in  Figure  3.1.  We  see  that  for  the  more  viscous  finger,  there  is  reen- 
forcement  of  the  basic  finger  growth  on  the  upstream  side  and  a  cancellation  on  the 
downstream  side.  For  the  solvent  finger  these  trends  are  reversed. 


H 


O 


H 


Figure  3.1.  The  one  dimensional  noninleracting  layer  ansatz.  Schematic  finger  shapes  and 
pressure  fields  are  displayed.  The  dash  lines  represent  the  unperturbed  pressure  fields 
and  the  solid  lines  rq)resent  perturbed  pressure  fields.  H  and  L  denote  the  high  and  low 
pressure  sides  of  the  finger  respectively. 


The  result  of  this  picture  is  that  the  oil  finger  is  transported  downstream,  and  grows 
longer  on  its  upstream  side  only.  The  solvent  finger  has  only  very  slow  transport,  and 
grows  on  its  down  stream  side  primarily.  Funhermore,  the  Taylor-Saffman  instability  of 
the  solvent  finger  contributes  to  its  down  stream  growth  while  the  Taylor-Saffman 
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stability  of  the  oil  finger  yields  a  very  reduce-d  tendency  for  down  stream  growth.  For 
this  reason,  we  get  narrow,  long  fingers  of  oil,  and  broad  fingers  of  solvent.  The  solvent 
fingers  rapidly  approach  an  equilibrium  configuration,  in  which  the  relative  width  of  the 
two  channels  adjusts  until  the  condition  MV  =  1  is  achieved.  The  oil  fingers  move  down 
stream,  allowing  the  instability  to  re-establish  itself  near  the  inlet,  and  new  fingers  to 
form,  in  a  quasi-periodic  or  chaotic  manner. 


t=100 


t  =  400 


Figure  3.2.  Excess  in  flow  of  oil  in  the  bottom  channel  causes  the  growth  of  the  oil  finger 
into  the  solvent  channel.  We  have  shown  two  time  steps  for  the  growing  finger  at  /  =  100 
and  t  =  400.  The  mesh  size  is  70  x  30  for  the  elliptic  (pressure)  equation. 
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t  =  400 


t  =  600 


Figure  3.3.  Excess  in  flow  of  solvent  in  the  bottom  channel  causes  the  growth  of  the  sol- 
vent finger  into  the  oil  channel.  We  have  shown  two  time  steps  t  =  400  and  t  =  600  for 
the  growing  finger.  The  mesh  size  is  41  x  17  for  the  hyperbolic  equation  and  70  x  30  for 
the  elliptic  equation. 


Next  we  examine  the  solurion  of  the  full  Hele-Shaw  equations  in  these  two  cases. 
In  Figure  3.2  we  show  time  frames  presenting  the  interface  positions  for  excess  oil 
inflow,  while  Figure  3.3  treats  excess  solvent  inflow.  We  see  a  pattern  of  interface 
development  as  predicted  above.  To  continue  this  analysis,  we  show  in  Figure  3.4  the 
pressure  fields,  plotted  as  a  function  of  x,  along  fixed  values  of  y,  at  constant  t.  These 
figures  should  be  compared  to  the  qualitative  discussion  given  at  the  beginning  of  this 
section.  Then  we  compare  the  y  component  of  the  velocity,  also  plotted  vs.  x,  along 
y  =  const,  lines  at  the  same  fixed  t  in  Figure  3.5. 


11- 


Figure  3.4.  We  plot  pressure  fields  vs.  x  along  two  horizontal  lines  at  dimensionless 
heights  (i)  ylH  =  0.2  and  (ii)  >  ///  =  0.3  above  the  /  =  0  interface  position,  (a)  Excess  oil 
inflow,  t  =  100.  Inside  the  finger,  the  slope  of  the  pressure  curve  is  steeper  than  outside, 
(b)  Excess  solvent  inflow,  i  =  600.  Inside  the  finger,  the  slope  of  the  fH^ssure  is  smaller 
than  on  the  outside.  Also  the  slope  of  the  pressure  is  smaller  in  front  of  the  finger,  indi- 
cating a  stagnant  region  in  front  of  the  finger.  The  small  slope  of  ciuve  (i)  near  x  =  0  and 
the  non  monotone  portion  of  (i)  in  Figure  3.4a  is  due  to  the  influence  of  the  small  barrier 
separating  the  flow  channels  near  x  =0. 


12 


Figure  3.5.  Plots  of  v^  vs.  x  along  horizontal  lines  and  times  coinciding  with  the  plots  of 
Figure  3.4.  The  discontinuity  in  the  plot  occurs  at  the  interface  between  the  oil  and  the 
solvent. 


The  results  we  have  shown  from  Figures  3.2  to  3.5  are  not  sensitive  to  mesh 
refinement.  We  obtained  almost  identical  results  after  doubling  the  number  of  the  grid 
lines. 
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4.  Pinch  Off  and  Bubble  Generation  in  a  Uniform  Channel 

In  general  terms,  we  find  that  for  MV  >  1,  (2.2)  gives  Yo  >  0.  The  system  deviates 
from  equilibrium,  and  the  excess  oil  in  flow  produces  a  transverse  velocity  field  which 
drives  the  oil  into  the  solvent  channel,  as  shown  in  Figure  2.2.  To  describe  the  physical 
process  of  pinch  off,  there  are  two  natural  time  scales  to  consider.  The  first  is  the  pinch- 
off  time  tp,  which  is  the  characteristic  time  for  a  single  finger  to  grow  to  the  height  of  the 
channel  width  and  the  second  is  the  time  t^,  which  characterizes  the  expansion  of  the 
finger.  To  study  tp,  we  approximate  the  shap>e  of  the  finger  by 

xo(t) 


,     „      4/1  (r) 


3n 


ix-xo{t)f+xl{i) 


x>0     ,  (4.1) 


where  A  (t)  is  the  area  of  the  finger.  By  conservation  of  mass  and  incompressibility,  A{t) 
is  equal  to  the  total  flux  of  the  nonuniform  flow  v  -  v'  =  yv'  into  the  channel  with  excess 
flow, 

Ait)  =  WHy^v/t    .  (4.2) 

Here  H  is  the  height  of  the  entire  flow  channel  and  xq  defines  a  half  width  for  the  finger. 
Since  the  area  covered  by  yix,  t)  in  (4.1)  in  the  interval  0  <  j:  <  j:o  is  one  third  of  total, 
we  define  j:o(0)  as  follows:  the  nonuniform  flux  v  -  v'  crossing  this  interval  at  r  =  0  is  the 
same  ratio.  Therefore  in  case  of  large  M  according  to  the  theory  of  §2,  we  have 

':^^\+cxp{-Kx/2WH)_^ 

I      \-c\pi-Kx/2WH)  1 


7^^1-hexp(-7u/2\y//)^ 


3 


txpi-Kx/2WH) 

By  this  definition  of  Xq,  we  find  that 

X(i(.0)  =  0A15WH     . 

Since  the  oil  finger  penetrates  into  the  rapid  moving  solvent  channel,  its  width  expands 
with  time.  As  an  approximation,  we  assume  that  the  expansion  rate  is  constant  with  time 
and  profKjrtional  to  the  uniform  solvent  velocity.  Therefore 

;co(r)  =  0.175W//-(-;v',f  ,      0<C<1,  (4.3) 

where  ^  is  a  dimensionless  parameter.  We  found  that  adjusting  ^  =  C,{W)  to  depend  on  W 
only,  C,  =  0.2,  it  is  possible  to  obtain  a  good  fit  for  simulations  covering  a  wide  range  of 
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values  of  V,  M  and  W.  The  height  of  the  finger  is  given  by 

,(,) = ^m. . 'HM^L (4.4) 

In  Figures  4.1,  4.2a  and  4.2b  we  compare  the  interface  and  height  of  the  oil  finger  calcu- 
lated by  (4.1)  and  (4.4)  with  the  results  of  numerical  simulations.  The  equation  (4.4) 
defines  a  threshold,  or  maximum,  for  the  finger  height,  hj,  which  is 

Pinch  off  occurs  only  when  the  solvent  channel  width  is  less  than  the  finger  threshold, 
see  Figure  4.1b.  Therefore  we  have  a  condition  for  occurrence  of  pinch  off, 

(.-..X^L.         .  ,4.6) 

The  pinch  off  time  tp  is  defined  by  setting  the  finger  height  (4.4)  to  the  height  of  the  sol- 
vent channel  /j^  =  (1  -W)H,  yielding 

Q.\15WHM{\-W) 


tp  =  -r ^ \ .  (4.7) 

in 

When  the  pinch-off  occurs  repeatedly,  it  turns  a  laminar  flow  of  solvent  into  slug 
flow.  To  estimate  approximately  the  slug  size,  we  introduce  the  dimensionless  length 
scale 

a=-^.  (4.8) 

Over  a  range  of  simulations,  we  find  this  length  scale  to  characterize  the  size  of  the  initial 
slug  within  50%. 
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Figure  4.1a.  The  variation  of  finger  height  with  time.  The  solid  line  and  dashed  line  are 
the  numerical  and  theoretical  results. 
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t  =tp 


Figure  4.1b.  The  threshold  finger  height  hr-  When  the  solvent  channel  height  h,  satisfies 
hs  <  hr,  pinch  off  occurs.  When  hs  >  hr  pinch  off  does  not  occur. 
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t  =  20 


t  =  40 


t  =  60 


Figure  4.2a.  The  interface  at  limes  t  =  20, 40, 60,  befwe  and  near  pinch  off,  are  calculat- 
ed by  numerical  simulation,  for  V  =  2,  A/  =  5  and  W  =  0.82.  The  mesh  sizes  are  70  x  40 
hyperbolic  and  150  x  70  elliptic. 
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t  =  20 


t  =  40 


t  =  60 


Figure  4.2b.  The  interface  calculated  by  the  theoretical  model  (4.1)  at  the  same  times  and 
parameter  values  as  were  used  in  the  simulation  of  Figure  4.2a. 


The  second  time  scale,  r^,  is  the  characteiisric  time  for  a  finger  moving  down  stream 
to  expand  its  width  by  an  amount  2WH  lit, 

2WH/n 


ts  = 


CV. 


(4.9) 


To  give  an  order  parameter  approximately  characterizing  the  flow  pattern,  we  introduce 
the  ratio  P  of  the  two  time  scales. 


P=-  = 


3ji 


-  ;(1  -  W)M 


0.27CA/(1-H') 


(4.10) 


For  P  >  1 ,  the  finger  pinches  off  before  it  moves  far  enough  downstream  for  new  fingers 
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to  form.  The  consecutive  pattern  of  pinch  offs  turns  the  flow  stream  into  a  series  of 
slugs.  For  p  <  1,  the  rate  of  finger  growth  on  the  interface  is  slowed  down  by  its  expan- 
sion. The  fluctuations  of  the  interface  become  increasingly  significant  and  pinch  off,  even 
it  occurs  eventually,  may  not  be  predicted  accurately  by  our  model.  When  P  <  0  the  sol- 
vent channel  width  is  above  the  threshold,  pinch  off,  if  it  occurs,  is  a  random  event  gen- 
erated by  the  fluctuations.  In  this  case  the  pinch  off  time  is  not  well  defined.  We  general- 
ize (4.10)  to  classify  the  resulting  chaotic  flows  by  defining  the  amplitude  5/i  of  the 
fluctuation  required  for  an  oil  finger  to  break  the  solvent  channel, 

5/2  =  (l-W)--^      . 

Then  (4.10)  can  be  rewritten  as 

P  =  -0.27(1-VV)=^-'^"-^^'     ^  =  ^=3kCA/(1-U^)       •  ^'-''^ 

In  the  limit  of  laminar  flow,  we  have  MV  =  I  and  0  =  yQ  =  hj=x  and  so  P  =  -  3.7.  In 
summary  we  conclude  that  the  laminar-slug  flow  transition  occurs  near  p  =  0  and  in  this 
region,  the  flow  pattern  is  chaotic,  with  random  pinch  offs  or  multi-finger  growth  on  the 
interface.  For  large  positive  P  well  developed  slugs  are  formed.  When  P  approaches 
-3.7,  laminar  flow  results.  Our  simulation  results  indicate  a  qualitative  agreement  with 
this  prediction  for  the  nature  of  the  phase  plane.  We  summarize  our  simulation  results  in 
Table  4.1.  In  the  laminar  and  chaotic  flow  limits,  ^  is  not  well  defined  by  our  fitting  pro- 
cedure, and  we  have  arbitrarily  set  ^  =  0.2. 

For  MV  <  1,  i.e.  Yo  <  0.  there  is  excess  solvent  inflow.  In  the  region  of  physical 
interest,  V  is  positive,  so  IYq  I  can  not  be  too  large.  Moreover,  the  analysis  of  §3  shows 
that  in  this  case,  the  fingers  are  transported  slowly,  grow  on  the  down  stream  side,  and 
the  overall  flow  pattern  converges  to  a  steady  state  equilibrium  of  laminar  flow.  As  result 
the  pinch-off  is  unlikely  in  this  case. 
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V 

M 

W 

P 

O-th 

'J^mtm 

V 

^'P 

Flow  pattern 

0.11 

10 

0.82 

-0.34 

- 

- 

- 

- 

Laminar  flow 

1 

5 

0.2 

-1.4 

- 

- 

- 

- 

Chaotic  flow 

2 

5 

0.4 

0.1 

- 

- 

- 

- 

Pinch  off 

1 

10 

0.92 

4 

0.25 

0.20 

20 

20 

Slug  flow 

0.33 

15 

0.92 

2.1 

0.64 

0.55 

31 

28 

Slug  flow 

1 

5 

0.92 

3 

0.31 

0.4 

39 

37 

Slug  flow 

1.43 

15 

0.82 

2 

0.24 

0.15 

35 

33 

Slug  flow 

2 

5 

0.82 

2 

0.23 

0.15 

54 

58 

Slug  flow 

Table  4.1.  Values  of  slug  size  and  flow  regime  order  parameters  a  and  3  are  tabulated 
for  various  simulations,  to  display  the  dependence  on  the  velocity,  mobility  and  width  ra- 
tios V^,  M  and  W.  o^,  and  a^^  are  the  theoretical  and  numerical  values  for  slug  size  in 
units  of  channel  flow  width.  The  theoretical  and  numerical  values  of  the  pinch  off  times, 
tp  and  t'p,  respectively,  are  compared.  For  W  =  0.82  we  choose  ^  =  0.21  and  for 
IV  =  0.92,  C  =  0.16. 


The  numerical  calculations  of  initial  finger  growth  in  the  slug  flow  regime  are 
not  sensitive  to  the  grid  size,  and  converge  under  mesh  refinement.  At  later  times 
in  the  simulation,  i.e.  for  the  formation  of  multiple,  repeated  slugs,  and  their 
movement  down  stream,  we  find  quahtatively  similar  pictures,  differing  in  detail, 
as  the  mesh  is  refined.  It  appears  that  the  use  of  an  unregularized,  unstable  formu- 
lation of  the  physics  is  the  cause  and  that  the  use  of  siuface  tension  or  other  regu- 
larization  would  be  required  to  obtain  convergence.  See  Figures  4.3  and  4.4. 
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(a) 


(b) 


Figure  4.3.  Consecutive  pinch  off  turns  a  laminar  flow  into  slug  flow.  Here 
V  =  1,  A/  =  5,  W  =0.92  and  P  =  2.  We  compare  distinct  mesh  sizes,  (a)  The 
mesh  size  is  42  x  14  for  the  hyperbolic  equation.  84  x  30  and  for  the  elliptic 
equation,  (b)  We  repeat  the  calculation  of  (a)  by  using  a  refined  mesh,  which  is 
64  X  28  for  the  hyperbolic  equation,  180  x  60  and  for  the  elliptic  equation. 


(a) 


(b) 


Figure  4.4.  Comparison  of  distinct  flow  regimes,  (a)  The  pinch  off  occurs  in  a 
random  way.  Here  V'  =  2,  M  =  5,  W  =  0.4  and  P  =  0.1.  (b)  Chaotic  flow  without 
pinch  off.  Here  V  =  \,  M  =5,W  =0.2  and  P  =  -1.4.  The  mesh  size  for  both 
cases  is  30  x  10  for  the  hyperbolic  equation,  60  x  30  and  for  the  eUiplic  equation. 
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5.  Variable  Channel  Width  and  Complex  Barriers 

As  we  show  in  Figure  5.1,  even  if  the  boundary  conditions  give  steady  lam- 
inar flow  i.e.  even  if  the  velocity  ratio  is  equal  to  its  equilibrium  value,  the  flow  in 
the  channel  throat  can  be  out  of  equilibrium.  In  this  case  an  instability  of  the  inter- 
face develops  locally.  Since  the  width  ratio  W  and  velocity  ratio  V  change  with 
variation  of  the  channel  width,  3  may  change  in  a  variable  channel.  In  a  region 
where  p  is  positive  and  not  too  large,  pinch  off  is  favored,  leading  to  a  transition  of 
laminar  to  slug  flow.  The  point  to  our  analysis  is  that  the  pinch  off  can  occur 
upstream  of,  i.e.  upon  entering,  the  channel  throat. 


:7 


n 


Figure  5.1.  Both  the  inlet  and  outlet  injection  rates  are  set  to  their  equilibrium  values,  i.e. 
MV  =  1,  where  M  =  5.  However,  two  barriers  in  the  middle  produce  a  narrow  flow 
channel  and  alter  the  velocity  ratio,  so  that  a  change  of  flow  pattern  results.  At  the 
upstream  side  of  barriers,  V  increases  from  0.2  to  1.0  and  value  of  p  is  1.1.  A  transition 
from  laminar  to  slug  flow  is  observed.  The  mesh  size  is  60  x  20  for  the  hyf)erbolic  equa- 
tion, and  120  x  60  for  the  elliptic  equation. 


In  the  problem  of  two  viscous  fluids  flowing  among  random  arrays  of  bar- 
riers, the  complex  flow  channels  can  give  rise  to  competing  flows  which  causes 
pinch-off  as  we  discussed  in  previous  sections.  An  idealized  geometry  which 
displays  this  effect  is  shown  in  Figure  5.2.  First  the  solvent  inflow  is  split  by  two 
front  barriers,  and  then  driven  into  a  narrow  channel.  Because  of  the  competing 
flow  channel  around  the  first  barrier,  the  velocity  ratio  at  the  entrance  to  the 
second  channel  deviates  from  equilibrium  values,  i.e.  ^fV  ^  1.  Thus  pinch  off 
results  and  slugs  are  produced.  In  Figure  5.3,  we  show  the  results  of  numerical  cal- 
culations of  flow  in  a  random  array  of  barriers,  perhaps  more  representative  of 
reality.  In  this  flow,  pinch  off  appears  to  be  associated  with  competing  flow  paths, 
as  in  Figure  5.2. 
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Jk 


Figure  5.2.  Numerical  calculation  of  viscous  flow  of  two  fluids  in  a  chan- 
nel with  complex  geometry.  The  computation  shows  formation  of  com- 
peting flows  and  pinch  off  between  barriers.  We  use  gray  scale  to  distin- 
guish three  different  phases.  The  barrier,  oil  and  solvent  are  shaded  by 
black,  gray  and  white. 
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Figure  5.3.  Miscible  flood  of  an  oil  core  containing  an  array  of  imperme- 
able barriers  approximating  the  microstructure  of  porous  media.  The 
boundary  between  the  injected  solvent  phase  and  the  displaced  oil  phase  is 
shown  at  three  different  dimensionless  times  during  the  flood,  from  left  to 
right  t  =  1.1,  3.0  and  4.0.  The  mobility  ratio  M  =  10.  The  competing  oil 
and  solvent  flows  cause  pinch  off.  The  barreiers,  oil  and  solvent  are 
shaded  with  black,  gray  and  white. 
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6.  Conclusions 

We  have  studied  the  viscous  flow  of  two  fluids  in  a  Hele-Shaw  cell  of  com- 
plex geometry.  We  have  shown  that  when  the  velocity  ratio  of  the  two  fluids  is  not 
equal  to  its  equilibrium  value  i.e.  AfV  ^\,  a  new  type  of  instability  (i.e.  fingers 
growing  which  have  qualitative  features  not  present  in  the  classical  Taylor- 
Saffman  instability)  develops.  Using  a  simple  theoretical  model,  we  have  studied 
the  pinch  off  phenomena  quantitatively.  We  are  able  to  give  the  condition  for  the 
occurrence  of  pinch  off  and  calculate  pinch  off  time.  An  order  parameter 
P(V,  M,  W)  is  given  to  characterize  the  distinct  flow  pattern  regimes  and  the  tran- 
sitions between  them.  We  obtained  agreement  between  analytic  and  simulation 
results.  The  transition  of  laminar  to  slug  flow  provides  a  possible  mechanism  for 
foam  generation  in  enhanced  oil  recovery. 


7.  Appendix.  The  transverse  velocity  field  of  the  nonuniform  flow 

To  calculate  the  transverse  velocity   field  of  the  nonuniform  flow  along  the 
interface  at  t  =  0,  we  stan  from  the  elliptic  equations  (A.l) 


VV,  =0  ,      v,=-—VP,  ,    J  =  1,2 


(A.l) 


with  the  boundary  conditions 

Vx(.x  =  0.  y)  =  - 


^1. 

V2. 


v^(x=L,y)=' 


y  >0 
y  <0 


v'l.         y  >  0 
v'2,        y  <  0 


(A.2) 


v/x,  y  =  -WH)  =  Vy{x,  y=H-  WH)  =  0     . 

The  interface  is  located  y  =0,  at  r  =  0.  For  /:  ;t  0  the  eigen-solutions  of  (A.l)  have 
the  form 


P\  =  AkCOs[k(y  -  (1  -  W)H)]  cosh[)t(L  -  x)]  , 
Pi  =  Bi,coslk(y'  +  WH)]  cosh[A:(L  -  x)]       . 


(A.3) 
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For  /t  =  0 

P,  =  C,x,  J  =  1,2 

where  Ak,  B*  and  C,  are  constants.   The  condition  of  continuity  of  pressure  and 
normal  velocity  at  the  interface  gives 

P'{(y=0)  =  P^2(y  =  0)    , 
-^-l-P,(y=0)  =  ^^P2(y  =  0)     .  (A.4) 

From  (A. 3)  and  (A.4),  the  eigen-frequences  of  *  =ik„  can  be  determined  from  the 
equation 

sink„H  =  ^^     ^'  sin[A:„(l  -  2W)H]    .  (A.5) 

1^2  +^1 


The  solution  of  (A.l)  under  the  boundary  condition  (A. 2)  can  be  written  as 

2(n,6vi-H25v2)  1    .  2„  „,„, 


cosh  [k..(L  -x)] 

X  sin[^„(l  -  W)H]cos[k„y  -  k„(\  -  W)H] r^— (A.6) 

sinh  KnL 


and 

^2  =-M2V2X  -5:^2 .    2.,    ',7,'"     '       rV/,\\       ..,.„.  ^S'"'[^n(l-W)//] 


2(H,5vi-^i25v2)  1      .2, 

Hisin2(^„W//)-t-^2Sin2[*n(l  -W)//]  kl 


cosh  [k„{L  -x)] 

X  sm(k„WH)cos[k„(y  +  WH)] .   "  ,   ,  (A.7) 

sinh  k.L 


where 

Sv,  =  V,  -  v',     . 

We  obtain  the  expression  for  transverse  velocity  Vj, 

2ai,5vi-^25v2)  1 


Vy(x,y=0)  =  -'£- 


„  \iism^(k„WH)  +  \i2sm^[k„(\  -  W)H]  k„H 


2  2/ cosh  [k„(,L-x)] 


xsm\k„WH)sm^(k„{l-W)H) — — .  (A.8) 

sir)hk„L 
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When  L^oo,  we  have 


v/x,  >=0)  =  -X 


2()a.i5vi  -(i25v2) 


^isin2(^„W//)  +  \i2sm^[k„i\  -  W)H]  k„H 


X  sm^(k„WH)sm^ik„{\  -  W)H)c\p  (-k„x) 


(A.9) 


We  discuss  the  solutions  of  (A.9)  in  some  special  cases. 
1 


i)  When  W  =  —,  from  (A.5)  we  have 


sink„H  =  0  ,        k„  = 


nn 
H 


s\nik„WH)  =  sin 


nn 


sin(k„(\-W)H)  =  sm 


nn 


0,  n  =  2k 

(-1)*,  n=2k  +  \ 

0,       n  =  2k 

(-1)*,  n=2k  +  l 


(A.IO) 


and  thus 


v,=- 


2(^1 5v,  -^l25v2) 


1 


^l+^l2       „(2n  +  i)7c 

-In 


exp  [-{2n  +  1)7U/,//] 


(^tl6vl-^25v2)^  l+exp(-7U///) 
7i(|ii  +  |i2)  1  -  exp  (-7U:  ///) 

ii)  For  p.2  ^  M-i .  from  (A.5)  we  have 

sinJk„//=-sin)k„(l-2W)// 


so  that 


k„  = 


{In  +  1)71 
2WH 


We  obtain 


2(ni5vi  -^25v2)  2^,        1  ,  n  ,       2W 
cos  "^ 

^ll  +^2Cos2[(«  +  y)-^] 


(A.ll) 


(A.12) 


-cos^[(«  +  —)— ] 


2    W    (2n  +  l):: 


exp  [-(2n  +  1)7C/2W//] 


28- 


^   n       \-exp{-Tix/2WH) 
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